FINITE ELEMENT MODELING OF TRUS3 STRUCTURES j f 
WITH FREQUENCY-DEPENDENT MATERIAL DAMPING I 


GEG.iGE A. LES1EUTRE 

Aerospace Engineering 
The Pennsylvania State University 


4th WORKSHOP ON COMPUTATIONAL CONTROL 
OF FLEXIBLE AEROSPACE SYSTEMS 

Williamsburg, Virginia 
July 13, 1990 


FINITE ELEMENT MODELING OF TRUSS STRUCTURES 
WITH FREQUENCY-DEPENDENT MATERIAL DAMPING 

ABSTRACT 

In advanced engineering systems such as large space structures (LSS) or robots, the combination of severe 

knowledge of structural dynamic behavior, including the magnitudes and mechanisms of in eren P S- 

Material damping is likely to be an important, perhaps dominant, contributor to 
structures and to on-orbit damping in precision spacecraft. The damping of most engine g 
a strong dependence on frequency. 

A physically-motivated modelling technique for structural dynamic analysis that accomodates frequency- 
dependent ma"™ 1 ! damping has been developed. Key features of, lh« ™ .»5.1te 

augmenting thermodynamic Fields (ATF) to interact with the usual mechanical t,e “* ! method 

treatment of the resulting coupled governing equations using finite element analysis methods. The ATF method 
is fully compatible with current structural finite element analysis techniques. 

The method is demonstrated in the dynamic analysis of a 10-bay phnar truss ‘ ^JommodTSyses of 
representative of those contemplated for use in future space systems. Ana y 1 , , » -rp analysis 

and an undamped ttuss structure agree well in terms of modal but 0* 

also yields modal damping and complex mode shapes. The undamped feqtmnc.e m lower by as 
percent in higher modes because the relaxed (stat.c) modulus value ts used as ,«1 praettem m a 
preserving the characteristic frequency dependence of matenal ttorrplltg, the ATF method reflects the treque y 
dependence of material modulus as well. 

With the continued development of better analytical tools such as this ATF ! e^e^cwssible to 

modelled more accurately in the design of engineering systems and may ultimately become more 

design specification. 
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matrix fiber 


MATERIAL DAMPING WILL BE IMPORTANT IN PRECISION SPACE SYSTEMS 


Hi h , . en f neenng s y sterns such as large space structures (LSS) or robots, the combination of severe 

SS&S? flL St P"® ent requirements, and structural design constraints can result in structures which exhibit 

, ™ flexlb ‘! lt >’' Passive and active damping of these structures is important for several reasons. In terms 

w h ‘ sher da , mpl f 8 u Can reduce stead y- state vibration levels and can reduce the time needed for 
transient vibrations to settle Inherent passive damping can reduce the magnitude of control needed, and can 

In f ^ SyStCm c . omplexit >' Passive dam P' n g can also strongly couple vibration modes which are closely- 
rTcnnnH frequ f n< ; y and C0I ? puted assuming no damping. Most importantly, however, the design of stable, fast- 

benefi ' s f ™ m ac ™ rate knowledge of structural dynamic behavior, which depends on 
tne magnitudes and mechanisms of inherent damping. F 

for sucS^H^r* da f mping exist / or a 8' ven aerospace vehicle. A potential classification scheme 
in sipnm^ni d shown above. In general, damping is not well-quantified in the design process, which results 
f” madC ° n thC baS ' S 0f qualit ‘ tative information. There is considerable room 
n P . a ! n e analysis of damping in aerospace systems. For example, a damping design budget might 
allocate contributions to total damping required from individual damping sources. 

A case can be made for the importance of material damping in at least some precision space vehicles. In 

svsTem° damning S f ( , mCtlJreS wh, ^ h operate in , the , atmosphere, air damping and joint damping typically dominate 
.r/h P 8 ' r ’ air damping 1S clearl y eliminated in space, and the effects of joint damping will be 
of requirements for precision ("tight" joints) and typically low vibration levels (friction 
n P. „ a ena dam P‘ n S 1S thus likely to be an important, perhaps dominant, contributor to damping in 
monolithic structures and to on-orbit damping in precision spacecraft. 
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w a erial damping is comply: 


• DEPENDS ON FREQUENCY, TEMPERATURE, DEFORMATION, 
AMPLITUDE, GEOMETRY 



. ATF METHOD ADDRESSES FREQUENCY-DEPENDENCE 


Material damping is generally a complex function of frequency, temperature, type of deformation, 
amplitude, and structural geometry. 

The figure on the left is adapted from the frontispiece of the pioneering text, M tt and AnelasticitY Jai 
Metals and illustrates the typical frequency-dependence of material damping. Note that there are a 
discrete damping peaks, and that a physical (possibly atomic) process is associated with each peak (e.g.. 
"transverse thermal currents”). 

The figure on the right shows some data recently obtained for aluminum in a flexural vibration test. The 
frequency dependence of damping is apparent. Similar data were obtained for graphite/aluminum composit 
materials, although weaker frequency dependence was observed. 

Current popular treatments of damping in structural dynamics are generally not physically-motivated and 
do not preserve the fundamental frequency-dependence of material damping. This observation provided the 
motivation for the subject work. 
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Df MPIN'i MODELS AND RELATED WORK 
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Several methods are currently used to incorporate the effects of material damping into structural models. 
These methods include viscous damping, frequency-dependent viscous damping, complex modulus, hysteretic 
damping, structural damping, viscoelasticity, hereditary integrals, and modal damping. Modal damping is 
probably the most widely used in structural control design applications. Each has some utility, but each suffers 
from one flaw or another. For example, a one-dimensional structure made from a single material with viscous 
damping would exhibit damping that increases monotonically with frequency — such behavior is not observed in 
practice. Although some damping models, such as viscoelasticity, have the potential for better accuracy than 
more widely-used methods, they are not commonly used in the engineering community — perhaps because of the 
lack of physical motivation, difficulty of use, or lack of data. 

Other researchers have noted the inadeqacy of current damping modeling techniques and have addressed 
the development of improved methods. 

Golla, Hughes, and McTavish (GHM) of the University of Toronto have developed a time-domain Finite 
element formulation of viscoelastic material damping. Their work was guided by the observation that 
experimental results, often recorded in the frequency domain, are of little direct use in time-domain models. 
Their results resemble those reported here in an important way — in the introduction of additional "dissipation 
coordinates.” It might not be surprising that additional degrees of freedom would be required to model material 
behavior with increased accuracy. However, the GHM approach is fundamentally a mathematical one, 
developing time-domain realizations from frequency-domain models — no attempt is made to provide a physical 
interpretation of the dissipation coordinates as thermodynamic field variables with a direct relationship to 
microstructural features of real materials. The GHM technique has been successfully used to fit a portion of an 
experimentally-determined curve of damping versus frequency, and standard analysis tools can be used to solve 
the resulting equations. 

Bagley and Torvik of the Air Force Institute of Technology have also developed a relevant model of 
material damping. The core of their concept is the use of fractional time derivatives in material constitutive 
equations. Their development was motivated by the observation that the frequency dependence observed in real 
materials is often weaker than the dependence predicted by first-order viscoelastic models. With four and five 
parameter models, they have been able to accurately represent the elastic and dissipative behavior of over 100 
materials over frequency ranges as broad as 8 decades. The application of the general fractional derivative 
approach to time-domain analysis is an area of continuing research. 
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E.ASIS FOR ATF 


. "INTERNAL STATE VARIABLES' OF MATERIALS SCIENCE MOTIVATE 
INTRODUCTION OF ' AUGMENTING THERMODYNAMIC FIELDS" (ATF) 


CONSTITUTIVE EQUATION 
ATF DYNAMICS 


a = £>- 5^ 


B 5 


• FREQUENCY-DEPENDENT CAMPING AND MODULUS 




(SMALL DAMPING) 


Structural dynamicists are the unintended beneficiaries of a sizable literature on material damping. For 
many years* crystallographers and metallurgists have used "internal friction" as a probe into the under ying 
structure of materials. By measuring damping as a function of frequency, temperature, load type, and 
amplitude, they can determine the mobility and activation energies of various microstructural features ot 
materials. These researchers have identified a multitude of internal variables and relaxation mechanisms which 
range in geometrical scale from crystal lattice dimensions on up to structural dimensions and, in temporal scale, 
over a similarly broad range. 

The "internal variables" of materials science suggest the introduction of augmenting thermodynamic fields 
(ATF) to interact with the usual displacement field of continuum structural dynamics. In the one-dimensiona 
case where stress is coupled to a single augmenting thermodynamic field (ATF), and where the ATF dynamics 
are governed by a first-order relaxation equation, the effective material modulus increases with frequency to an 
asymptotic value of E u , and a damping peak is observed. This is in accord with the observed be avior o many 
engineering materials. 

The ATF modeling approach essentially implements a continuum version of the "standard ^elastic sol “*‘ 
Weaker frequency dependence, such as that observed in high damping viscoelastic polymers, can be a esse 
through the introduction of multiple ATFs. 
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STRAW MAN PROELEM 


• MODAL ANALYSIS OF ATF-DAMPED TRUSS 

• ELEMENTS WITH DIFFERENT ELASTIC AND DISSIPATIVE PROP'S 


Longerons Diagonals Battens 



30 meters 


• COMPARE TO RESULTS OF CONVENTIONAL UNDAMPED ANALYSIS 


The ATF modeling approach is illustrated through application to the dynamic analysis of a large planar 
space truss, shown above. Such a structure resembles those proposed for many future space missions. It is an 
efficient beam-type structure built from 3 basic structural elements: longerons, which are parallel to the beam 
axis; diagonals, which bisect each rectangular bay; and battens, which are oriented transverse to the beam axis. 
The total length of the structure is 30 meters and the truss depth is 2 meters. 

Members of the truss of interest are connected through ffictionless pins so that no bending moments are 
transmitted through the joints, and so that no damping is introduced by the joints. Joint mass and other non- 
structural masses are ignored. Finally, each kind of structural element is assumed to be made of a diiferent 
material having different elastic, inertial, and dissipative properties. For this problem, it is assumed that a single 
augmenting thermodynamic field is sufficient to characterize the dissipative properties of each material in the 
frequency range of interest. 

C /wct nti0naI dam P in S analysis, if damping were considered at all, would likely employ the modal strain 
energy (MSE) method to estimate modal damping," This is an iterative process, requiring analysis of the 
undamped structure to determine the frequencies and mode shapes of the undamped modes. The mode shapes are 
used to determine the distribution of strain energy over the structure, while the frequencies are used to 
determine frequency-dependent elastic and dissipative properties. Modal damping is then estimated for each 
mode and, roughly speaking, is numerically equal to the sum of material damping ratios weighted by the fraction 
of the strain energy stored in each material. The MSE method can lead to serious errors when modes are closely 
spaced in frequency and when the damped mode shapes are much different than the undamped mode shapes. 

The vibration modes and frequencies of this ATF-damped truss are sought for comparison to those of the 
undamped case. A one-dimensional ATF-damped finite element appropriate for modeling the dynamic behavior 
of a single truss strut is developed first, then used in the analysis of the overall structure. 
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FINITE-LENGTH ROD 
EXECUTING LONGITUDINAL VIBRATION 



ENDS MECHANICALLY UNCONSTRAINED 


Consider the case of one-dimensional motion, corresponding to longitudinal vibration of a thin rod of 
length L. The mechanical displacement along the rod is denoted by u(x) (strain e(x)=u'(x)), and the rod has 
uniform mass density p and unrelaxed or dynamic modulus of elasticity E. A single augmenting thermodynamic 
field, i;(x), is introduced. 


The material property 5 describes the strength of the coupling of the two dependent fields, u and t In the 
absence of coupling of the two fields, increments of stress and strain are proportional, with E the relating factor. 
Analogously, a is the material property that relates changes in A, the thermodynamic conjugate of C, to those in 
5. B is the inverse of the time constant for uncoupled relaxation of the augmenting field. 
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GOVERNING EQUATIONS 


PDEs 


PRIMARY FORM 


pu~Eu" = 


f Bb ' 


a 




BCs 


ALTERNATE FORM 


pii-En" = -5y 
y+5y=| 


'35' 


a 

v j 


U " 


At ends (x=0 andx = L) 

u(x) = u 0 or a(x)=a 0 
ATF are internal variables (no BC) 


The equation of evolution for the mechanical displacement field, u(x), is developed from consideration of 
momentum balance (zero body forces are assumed). The equation of evolution for the augmenting 
^ ern u 0C ^ nam * C ^ ie ^’ ^ * S f° unc * t ^ rou §h th e use of a basic assumption of irreversible thermodynamics namely 
that the rate of change of £ is proportional to its deviation from an equilibrium value. This results in a first- 
order differential equation, a '’relaxation" equation. The result is a set of two coupled partial differential 
equations in u and £. 

The augmenting thermodynamic Field is essentially an internal field, Le. t there are no explicit boundary 
conditions that it alone must satisfy. However, the mechanical displacement field must satisfy either displacement 
( geometric ) or stress ("natural") boundary conditions at each end of the rod, as is the case in undamped 
structural dynamics. Note that the stress boundary condition involves the augmenting field, 

An alternative formulation of this one-dimensional case may be considered. For example, the preceding 
primary equations can be expressed in terms of y, the gradient of the £-field. Such a formulation contains only 
even spatial derivatives, and leads to some benefits in numerical solution, such as symmetric element submatrices 
and better convergence in terms of damping versus frequency. 
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FINITE ELEMENT TREATMENT 



u(x, t) = (. T (x) D 


U,(t) 

U 2 (t) 


7(x,t) =n T (x) T 


Y-(t) 

Y 2 (t) 


METHOD OF WEIGHTED RESIDUALS 


SINGLE ELEMENT 
MATRIX NOTATION 



ALTERNATE FORMS POSSIBLE 


M{i* + K{i 4 + F{it = 0 


The method of weighted residuals (MWR) is used to develop element matrices. The u-y formulation of the 
equations has been found to be superior to the u-£ formulation, and is shown above. Integration by parts is 
employed, changing the continuity required of the approximating and weighting functions. 

The same functions used to approximate the behavior of the dependent fields in the spatial region bounded 
by the element are used as weighting functions — when there is only one dependent field, this is known as 
Galerkin's method. 

The resulting sets of equations may be written compactly in matrix form, as shown. 

In this treatment, both dependent fields are approximated with linear interpolation functions. The 
corresponding element and the nodal values for the two dependent fields, u and y are shown in the figure. 
Anticipating solution of a first-order matrix eigenvalue problem, and to facilitate global matrix assembly, the 
elemental degrees of freedom are ordered as follows. 
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MATRIX EOUATIONS 



[ M]{ u} + [ K] { u } = { 0 } 

““lpl+f 0 , ^1(“L(°1 FIRST-ORDER FORM 

L 0 I J1 u] L - I 0 u / 1 0 J 



The structure of the undamped matrix equations of motion are shown in both second-order and first-order 
form. The structure of the ATF-damped matrix equations is shown for comparision. The ATF equations are 
more complex, describing the dynamics of the augmenting field as well as the coupling of the two fields. 

The ’’augmented mass matrix * 1 is generally symmetric and positive definite, while the "augmented stiffness 
matrix" is neither. Certain submatrices, K and H in particular, are generally symmetric and positive 
semidefinite, while the properties of the coupling submatrices, B anf F, depend on the specifics of the numerical 
approximation scheme employed. 
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ELEMENT MATRICES 
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For future reference and comparison, details of the elemental mass and stiffness matrices for the 
undamped and ATF-damped cases are shown. The matrices for the undamped case are shown in sjond-order 
form, while the ATF-damped matrices are shown in first-order form. Note that the elements of the undamped 
matrices also appear in the damped element matrices. 

Note that the material modulus, E, is different in the two cases shown. In the undamped case it is the 
static or "relaxed" modulus normally used in structural analysis. In the damped case it is the unrelaxed or 
high-frequency asymptotic modulus. The unrelaxed modulus is always greater than the relaxed value. 
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FREE V 5 3RATION PR^BLE^ 


• FREE LONGITUDINAL VIBRATION OF ROD 
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* GENERAL, UNSYMMETRIC EIGENVALUE PROBLEM 

* COMPLEX MODES 


In order to evaluate the performance of this formulation of an ATF-damped rod element, a specific 
problem was addressed, namely the determination of the natural modes of vibration of a free-free rod. 
Accordingly, no geometric boundary conditions were enforced. 

The elastic properties used correspond roughly to those of aluminum in SI units. The ATF and coupling 
properties were chosen to yield a peak damping ratio of 0.01 at the frequency of the 5 th mode. Note that 
numerical values for a and 5 cannot be uniquely specified in this approach. 

The matrix equations of motion were formulated and a general unsymmetric eigenvalue problem solved to 
yield complex eigenvalues, X, and mode shapes, x. Global system matrices are assembled from element matrices 
and geometric boundary conditions are enforced in the usual manner of structural finite element analysis. The 
damping ratio, £, for each mode was determined as the ratio of the negative of the real part of the eigenvalue to 
the total magnitude and plotted against the magnitude. 

Note that the spectrum of eigenvalues contains "vibration modes," "relaxation modes," and "rigid-body 
modes. 11 In the complex plane, the damped vibration modes lay near the imaginary axis, slightly in the LHP with 
negative real parts; the relaxation modes lie on the negative real axis. These relaxation modes are characteristic 
of the response of the y field. 
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RESULTS OF TRUSS ELEMcNT MCDA . A! ALYS'S 



u- Y EQ'NS PROVIDE 
BETTER DAMPING 
CONVERGENCE 


Typical numerical results for modal damping versus frequency obtained using this approach are shown 
above. The solid line shows the results expected on the basis of approximate Fourier analysis of the governing 
equations, while the symbols show the results of ATF-damped modal analysis using 10 elements. The square 
symbols indicate results obtained using the primary form of the governing equations, while the triangular 
symbols show the results of using the preferred, alternate form of the equations. 

Preservation of the characteristic variation of material damping with frequency is apparent in the results. 
As previously noted, conventional damping modeling techniques are incapable of producing such results. The 
frequencies predicted using this method appear to converge from above, as is the case with undamped elements. 
The finite element implemention based on the alternate form of the governing equations exhibits the proper 
variation of damping with frequency at all frequencies, while the implementation based on the primary equations 
Tolls off too rapidly at high frequencies. Predicted frequencies and damping from both implementations 
improve with increasing numbers of elements. 
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large: space truss modal analysis 



30 meters 



LONGERON 

DIAGONAL 

BATTEN 

CS area (m2) 

31 e-5 

19 e-5 

6.3 e-5 

Modulus (unrelaxed) (Pa) 

36.72 elO 

18.72 elO 

8.4 elO 

Density (kg/m3) 

2200 

1600 

2700 

Peak damping ratio 

0.005 

0.01 

0.05 

Frequency of peak (r/s) 

200 

2000 

8000 


The utility of this ATF modeling method is demonstrated in a modal analysis of the previously-described 
strawman 10-bay, 30-meter planar truss structure. The damped rod elements just described are modified to 
include the kinetic energy of transverse motion, and to address requirements for interelement continuity of the 
augmenting fields. 

The table above summarizes the key properties of the truss member materials. In particular, the peak 
damping ratios for each material are different, and are found at different frequencies. 
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RESULTS OF TRUSS MODAL ANALYSIS 






SOME 

UNDAMPED 
MODE SHAPES 


The figure above shows the results of the truss modal analysis m terms of damping versus frequency. 
Where conventional analysis using the MSE method would have required considerable effort to generate multiple 
modal damping values, the ATF method delivers them in a single, standard modal analysis. In addition, it 
delivers more accurate complex modes. 

Note that the damping in the lowest bending modes increases quite rapidly with frequency. This is due to 
two unrelated factors: First, in this kind of truss structure, the fraction of strain energy stored in the diagonal 
members increases with mode number-this is analagous to the effects of transverse shear in an isoropcbeam. 
Second, in this case, the damping of the material from which the diagonal members are made is increasing a 
frequencies below 2000 rad/sec. Also note that the 6* flexible mode is an extensional mode, and has 
'onsiderably lower damping than neighboring bending modes. 
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RESULTS — T RUSS MCD/.L. ANALYSIS 


MODE NO. 

' UNDAMPED 
FREQUENCY 

(R/S) 

DAMPED 

FREQUENCY 

(R/S) 

\ 

DAMPING 

RATIO 

(10e-3) 

4 

228.G 

231.4 

0.47 

5 

472.2 

480.8 

0.46 

' 6 

701.6 

719.3 

0.77 

7 

887.2 

917.5 

1.32 

8 

1024. 

1069. 

2.04 

9 

1102. 

1127. 

0.70 

10 

1109. 

1168. 

2.77 

J 


resuIt S F0 The AT^fnH a T da 7 ed a , naIysis was P erformed ' The table summarizes the frequency and damping 

first mode The undamnp/f” 156 reSU tS a ^ ee wed in terms of frequency, differing by about one percent for the 

er by as much as flve P ercent in the tenth mode because the 
relaxed (static) modulus value was used for that analysis, as is usual practice. 
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SUMMARY 


ATF METHOD 

\ ■ 


AUGMENTING THERMODYNAMIC FIELDS 


. PHYSICALLY-MOTIVATED, PRESERVES CHARACTERISTIC 
VARIATION OF M'.TERIAL DAMPING WITH FREQUENCY 


• FULLY COMPATIBLE WITH CURRENT FE TECHNIQUES 

• MIX DAMPED, UNDAMPED ELEMENTS 
. MULTIPLE FIELDS TO FIT DATA 


f FUTURE WORK ] 

. DAMPING BUDGET 
. CONTROLS 

• MSE COMPARISON 

• OTHER ELEMENTS 


. OPTIMIZATION 
. IDENTIFICATION 
. MATH ASPECTS 


* physically-motivated material damping model 

analysis methods has been developed. Termed the Augmeni g displacement field of continuum 

large space truss structure. 

This method preserves the characteristic method readily 

single augmenting field per material are readily extended P increase in the accuracy of a 

coordinates are required t0 represent 

additional aspects of material behavior, viz, damping, 
tmated™.*^ 
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